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ABSTRACT 

We present results of a comprehensive analysis of the X-ray pulsar X Persei over the 
period 1996 to 2011, encompassing the quite low state and subsequent strong outburst 
activity. Using data from the RXTE and Swift observatories we detected several con- 
secutive outbursts, during which the source luminosity increased by a factor of ~ 5 up 
to Lx — 1.2 x 10 35 erg s" 1 . Previously the source was observed in a high state only 
once. The source spectrum in a standard energy band (4 — 25 kcV) is independent of 
the flux change and can be described by a model that includes both thermal and non- 
thermal components. With the help of the INTEGRAL observatory data we registered 
the highly significant cyclotron absorption line in the source spectrum and, for the first 
time, significantly detected a hard X-ray emission from the pulsar up to ~ 160 kcV. 
We also report drastic changes of the pulse period during the outburst activity: a long 
episode of spin-down was changed to spin-up with a rate of P/P ~ — (3 — 5) x 10~ 4 
yr , that is several times higher than previous rates of spin- up and spin-down. To 
search for a correlation between the X-ray and optical lightcurves we took data from 
the AAVSO International Database. No significant correlation between optical and 
X-ray fluxes at any time lag from dozens of days to years was found. 
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1 INTRODUCTION 

4U 0352+309/X Persei is a classical persistent Be/X-ray bi- 
nary system, consisting of an X-ray pulsar and a Be star 
companion optically identified with the star HD 24534. It 
was discovered during a high X-ray intensity state in 1972 
(for simplicity below such high states are referred to as out- 
bursts), when pulsations with a period of ~ 835 s were de - 
tected by the Copernicus observatory (I White et al.lfl97^ ). 
The distance to the source was estimated by different au- 
thors in the range of 700 to 1300 pc. In t his paper we 
used t he value of 950 ± 200 pc obtained by iTelting et alj 
(|l998h . Adopting this distance the source peak luminosity 
Lx — 2 x 10 35 erg s _1 in the 2—10 keV energy range 
was registered in 1975. The close proximity allowed the 
source to be observed in X-rays by different missions and 
experiments after 1978, despite its relative low luminosity 
L x ^ (2 - 4) x 10 34 erg s _1 . 



E-mail: aal@iki.rssi.ru 



For a long time X Persei was th e only known persis - 
tent X-ray pulsar in Be-systems, until iReig fc Rochel (|l999h 
pointed out the existence of a subclass of binaries with Be 
companions which harbor a slowly rotating neutron star 
and are characterized by a persistent low-luminosity X-ray 
emission and, proba b ly, a long orbital period. Soon after, 
iDelgado-Marti et alj <|200lh succeeded in determining or- 
bital parameters for X Persei and showed that the pulsar is 
in a moderately eccentric orbit (e = 0.11) with a P or b — 250 
days orbital period. Recently, iTsygankov et all (|201 lh re- 
ported a similarly long orbital period of ~ 155 days for an- 
other X-ray pulsar RXJ0440. 9+4431, which belongs to this 
subclass as well. 

Usually spectra of X-ray pulsars have an exponential 
cutoff at energies about ^20 keV (see, e.g., the recent re- 
view of iFilippova et aill2005l . and references therein). How- 
ever results of the Ginga observatory gave evidence for the 
presence o f a hard tail in th e spectrum of X Persei above 
~ 15 keV (|Robba et al.lll99rj). Based on data of the Bep- 
poSAX observatory, Di Salvo et alj (|l998l ) showed that its 
broadband spectrum can be described by a combination of 
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the standard model and a power law with both high- and 
low-energy cut-offs, which dominat es above ~ 15 keV. Later 
lLa Palombara fc Mereehettil (2007) suggested that the ther- 
mal component with temperature of kT > 1 keV is typical 
for spectra of low luminosity, long period Be/NS binaries. 

Similar to other accreting pulsars, the history of the 
XPersei pulse period shows episodes of both spin-up and 
spin-down. Previous to 1978 its spin history was erratic, with 
a general spin- up trend at a rate oiP/P ~ — 1.5 x 10 -4 yr _1 . 
Since 1978 the neutron star has been s pinning down with a 
simil ar rate P/P ~ 1.3 x 10~ 4 yr -1 (fPelgado-Marti et all 
l200ll ). Note that the spin- up episode in 1975-1978 coincided 
with the decay phase of the outburst. 

Be stars are known to have a dense slow disk-like equa- 
torial wind along with a low-density fast wind at higher 
latitudes. XPersei has shown large variations in optical and 
infrared brightness (V = 6.1 — 6.8, K = 5.2 — 6. 7) and in the 
strength of the emission Balmer line (see, e .g. iRoche et al.l 
1 19931 ; belting et alj Il998l : IClark et al] |200ll . and references 
therein). These variations are usually interpreted as changes 
in the physical properties (density, geometry) of the Be star 
equatorial disk . The bright state of the star correlates with 
the presence of the emission lines and is evidence of the equa- 
tor ial disk presence. At the same time, as it was mentioned 
by ITelting et al.l (|l998T ). the strong variability seen in the 
optical and infrared is not reflected in the X-ray light curve. 
Particularly, the 1975 peak of the X-ray flux happened dur- 
ing the diskless state, while during the 1990 diskless state 
no X-ray brightening was observed. The absence of a di- 
rect correlation between the optical and X-ray light curves 
excludes the possibility that the neutron star is always ac- 
creting fro m the Be disc, pro bably due to the truncation 
of the disc (|Clark et al.ll200ll ). On the other hand the de- 
tailed study of the correlation between optical and X-rays 
was hampered by the absence of continuous measurements 
of the source flux in X-rays. 

In this paper we present a comprehensive analysis of the 
properties of the X-ray pulsar XPersei over the period 1996 
to 2011, encompassing the quite low state and subsequent 
strong outbursts activity. Thanks to the RXTE and Swift 
observatories, which have regularly monitored the whole sky 
in the last decades, we are able to look for the first time 
for a correlation between X-ray and optical data on a time 
scale from days to years. Moreover, based on data from the 
INTEGRAL and RXTE observatories we traced the evolu- 
tion of the source pulse period and its spectral parameters 
during outbursts. We also used deep observations of the IN- 
TEGRAL observatory to detect for the first time a hard 
X-ray emission from the pulsar up to ~ 160 keV at a high 
significance level. 



2 OBSERVATIONS AND DATA ANALYSIS 

In this work we use data obtaine d with the RXTE 
jBradt et al.lll993l ) and INTEGRAL jWinkler et all 120031 ) 
observatories. In particular, we used data of the 
ASM/RXTE monitor to trace the long-term behavior of the 
XPersei flux in the soft energy band (< 12 keV). To com- 
pare it with the source behavior in hard X-rays we used 
observa tions of the IBIS/IN TEGRAL and BAT /Swift tele- 
scopes (|Gehrels et alj|2004h . In the standard X-ray energy 



band (4 - 25 keV) data from the PCA/RXTE instrument 
(Obs. ID. 30099, 40424, 50404, 60067, 60068) were used to 
study the variability of the source spectrum and pulse pe- 
riod. 

In 2003-2006 XPersei was in the field of view of tele- 
scopes of the INTEGRAL observatory only episodically. 
During deep observations of sky fields around the source in 
2008-2010 the effective exposure of both X-ray telescopes 
of the observatory, JEM-X and IBIS/ISGRI, rose up to 
~ 147 ksec and ~ 3.7 Msec, respectively. The reduction of 
the ISGRI detector dat a was done usin g recen tly developed 
methods described by iKrivonos et al.l (|2010l ). The source 
spectrum from the JEM-X telescope was extracted with 
the standard 05^4 package version 9.o3- For timing anal- 
ysis the IBIS data was processed with the software devel- 
oped and maintained in the National Astrophysical Institute 
in Palermo^; a description of the d ata p rocessing technique 
can be found in lMineo et al.l (|2006l ) and lSegreto fc Ferrignd 
l|2007l ). 

To build a complete picture of the secular variability 
of the XPersei pulse period we included in our analysis 
data from the XMM-Newton observatory that observed the 
source twice in Feb 2003 and Feb 2010 (Obs. ID. 0151380101 
and 0600980101, respectively). Version 11.0 of the XMM- 
Newton Science Analysis System SAS was used to process 
the data. 

The final spectral and timing analysis for all X-ray 
instruments was done using the standard tools of the 
FTOOLS /LHEASOFT 6.7 package. 

To search for a correlation between the source outbursts 
activity in X-rays and optics (V-band), we took data from 
the AAVSO International DatabastQ. 



3 RESULTS 

3.1 X-ray light curves 

The XPersei light curve obtained with the ASM/RXTE 
monitor in 1996-2011 indicates that during the first years 
of observations the source intensity was low and relatively 
stable at the level of ~ 10 mCrab, that corresponds to a 
luminosity of L x ^ 2.3 x 10 34 erg s" 1 in the 2-10 keV 
energy band for the adopted distance (Fig[TJ). Such a level 
of intensity is typical for the normal (low) state of the pul- 
sar, which had been o b served since the en d of 1970s (see, 
e.g. IRoche et all 1 19931 ; ITelting et al. 1 1 19981 . and references 
therein). In 2000-2001 the source flux began to steadily in- 
crease and reached a maximum (~ 50 mCrab) in May 2003. 
It then dropped very quickly (within 4-5 months) by a factor 
of two and stayed near this level until 2008, when it began to 
increase again. The source flux at the new maximum (May 
2010) was approximately the same as the preceding one. And 
again it dropped very quickly by a factor of 2 (FigfTJ). 

The described source behavior at soft energies , as regis- 
tered by the INTEGRAL (20-60 keV) and Swift (15-50 keV) 
observatories (FigfTJ), almost completely mirrors the source 
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2 http://www.pa.iasf.cnr.it/ ferrigno/INTEGRALsoftwarc.html) 

3 http://www.aavso.org 
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Figure 1. (top) Light curves of X Persei obtained with the 
ASM/RXTE monitor in 1.3 - 12 keV (filled circles) and 
IBIS/INTEGRAL telescope in 20 - 60 kcV (open circles) en- 
ergy bands. The grey box indicates an interval with simultane- 
ous BAT/Swift observations, (bottom) Corresponding light curve 
of the source obtained with the BAT/Swift telescope in the 
15-50 keV energy band. To optimize the signal-to-noise ratio 
light curves were averaged over 30 and 20 days for ASM/RXTE 
and BAT/Swift, respectively, while ISGRI/INTEGRAL measure- 
ments were averaged over one spacecraft revolution (about three 
days). 
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Figure 2. Secular variabilit y of th e X P ers ei pulse period 
as his torically obse r ved by |Nagasel 1 19891); i Lutovin ov et al.l 



jl994h;lHaberiijl994 l ;lRobba et al.l il996T ); lDi Salvo et al.l J1998T) ; 
iDelgado-Marti et alj 1120011 1 fdark circles). Results of RXTE, IN- 
TEGRAL and XMM-Newton measurements from this paper are 
shown by the triangles, open circles and dark squares, respec- 
tively. 



variability in hard X-rays. The source outburst activity con- 
tinued during 2011, when a third outburst with a slightly 
lower intensity was detected by the Swift observatory. A 
comparison of X-ray light curves of X Perse i from Fig[T| 
with t h e source long-term b ehavior reported bv lRoche et al.l 
l| 1993ft l imiting et al. I (|l99Sh shows that a similar outburst 
activity was observed more than 35 years ago, in the middle 
of 1970s, with the Ariel V and Copernicus observatories. 
Moreover, the maximal flux from the source in 1975, was 
even slightly higher that measured in 2003 and 2010 and 
corr esponded to a sour ce luminosity of Lx ~ 2 x 10 35 erg 
s" 1 (|Telting et al.lll99&t ). 

Besides these large outbursts with durations of years, 
the light curves of X Persei demonstrate a tentative variabil- 
ity with the amplitude of about 5 — 10 mCrab. To search 
for a possible correlation of these source intensity varia- 
tions with the orbital phase, we folded ASM/RXTE and 
BAT/Swift light curves wi th a period of P or b = 250.3 days 
|Delgado-Marti et al.|[200ll ). but we did not find any signif- 
icant deviations from the constant either in soft or in hard 
X-rays. Finally, note that the intensity rise, which has been 
very smoothly for several hundred days for all observed out- 
bursts from X Persei, is non-typical for outbursts from usual 
transient sources (like black holes or transient X-ray pul- 
sars). It can be connected with other mechanisms or pro- 



cesses responsible for the outburst activity in the source un- 
der investigation. 

3.2 Pulse period evolution 

The secular variability of the X Persei pulse period is pre- 
sented in Fig[21 It includes results of measurements taken 
from the literature and calculations performed in this work. 
The photon arrival times have been preliminarily corrected 
to the barycenter of the Solar System. After that the epoch 
folding method has been applied to determine the period of 
the pulsar. The errors in th e period were es t imate d by the 
method briefly described bv lFilippova et al] (|2004l ). To cal- 
culate them we performed an analysis of 10 3 simulated light 
curves. Data in each light curve were generated arbitrarily 
inside the error bars of the original data points. After that 
the best period for each simulated light curve was deter- 
mined. Finally we computed the standard deviation of the 
best periods distribution and used this value as an estima- 
tion of the pulse period error. 

Due to the faintness of the source and a relatively long 
pulse period (~ 14 min) observations with an exposure of 
several hours (at least) are needed to determine the pulse 
period somewhat accurately. This was not a problem for 
INTEGRAL and XMM-Newton observations, where expo- 
sures were about several dozen ksec, but the typical dura- 
tion of one RXTE observation is only several ksec. There- 
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fore we grouped closely spaced RXTE observations into 
one set (when possible) and calculated the pulse period 
for each of them. Data of the XMM-Newton observatory 
were highly affected by event pile-up due to a very high 
count rate . Therefore in our analysis w e used the same ap- 
proach as lLa Palombara fc Mereghettil (120071 ) and consid- 
ered only data from the annulus area, excluding the cen- 
tral part of the source PSF. Applying the standard process- 
ing described in the SAS v. 11 cookbook we determined the 
pulse period of XPersei: 838.52 ± 0.20 and 835.94 ± 0.05 
sec for February 2003 and February 2010, respectively. The 
quality of the data during the last observation was much 
better than during the first one as its duration was longer 
(~ 31 ksec vs ~ 126 ksec), therefore the corresponding un- 
certainty is lower. The obtained value of the pulse period in 
February 2003 is slightly lower tha n the value reported by 
lLa Palombara fc Mereghettil (|2007m . There are two possible 
reasons for this difference: we used a more recent version of 
SAS, and we used the data from MPS camera s only. As men- 
tioned bv lLa Palombara fc Mereghettil (|2007m data of the pn 
camera were affected not only by photon pile-up, but dur- 
ing that observation the source was imaged close to a CCD 
gap. Finally, note that our value of the XPersei pulse period 
for February 2003 agrees well with the measurements from 
the INTEGRAL observatory, performed in April 2003 (see 
Fig®. 

Measurements from the INTEGRAL and XMM-Newton 
observatories indicate that sometime before the 2003 source 
flux maximum a deceleration of the neutron star rota- 
tion, which had lasted from 1978, changed to an accel- 
eration (Fig(2J. This repeats approximately the pulse pe- 
riod behavior during the outburst in 1974-1978, when the 
spin-up episode coincided with its declining phase, but 
started ~ 400 — 5 00 days after the luminosity maximum 
(IRoche et al. I ll993h . The average spin-up rate since 2003 
P/P ~ —3.6 x 10~ 4 yr _1 is higher than the one observed 
between 1974-1978. Moreover, measurements of the INTE- 
GRAL observatory suggest the presence of something like 
a plateau in the pulse period history during the time inter- 
val MJD 53000 — 54500, when the source flux was more or 
less stable between the first and second outbursts. In such a 
case the pulsar is spinning-up only during outbursts with the 
even higher rate of P/P ~ — 5x 10~ 4 yr _1 . The high spin-up 
rate can be probably connected with additional matter and 
angular momentum transfer to the neutron star due to the 
continuing source high state after the first outburst, and sub- 
sequent outbursts. This picture can be naturally explained 
in terms of the current theory of quasi-spherical acc retion in 
low-luminosity X-ray pulsars (|Shakura et al.ll2012h . In par- 
ticular, the observed spin-up/spin-down transitions corre- 
late with changes of the source luminosity, and their numer- 
ical values reflect the ratio of X-ray fluxes during spin-up 
and spin-down (see Fig.l and eqs. 56-58 in IShakura et all 
2012). Moreover, the plateau at the pulse period history is 
apparently connected with the equilibrium state of the pul- 
sar. Using the equation 
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Figure 3. Variability of spectral parameters and the 2 — 20 keV 
model flux of X Pcrsci during the rising part of the first outburst 
as measured with the PCA/RXTE instrument (see text for de- 
tails). 



Mie is the accretion rate in units 10 g/s, deriving from 
the source luminosity, and v& is the relative wind velocity 
(in units 10 s cm/c), we can estimate «g as ~ 200 km/s, 
which agrees with measured or suggested val ues (see, e.g., 
IWaters et "aT]|l988 | ; iDelgado-Marti et alj|200ll ). 

As noted by Lutovinov et al. ( 20041 ') . pulsations from 
the source are detected up to 100 keV. The pulse fraction 
increases from ~ 30% in the 2 — 10 keV energy band to 



where P eq — 837 s is the equilibrium period, ^130 is the neu- 
tron star magnetic moment in units 10 30 G cm 3 , which can 
be calculated from the measured magnetic field (see below) , 



~ 60% in the 20 — 60 keV energy ba nd, typical for X-ray 
pulsars (|Lutovinov fc Tsygankovll2009l ). 



3.3 Evolution of the spectral parameters 

Observations of XPersei with the PC A instrument of the 
RXTE observatory were performed several dozens of times 
in 2000-2003. They cover the low state and the rising phase 
of the first outburst very well. It allowed us to study a depen- 
dence of the source spectral shape variability on its intensity. 

Using RXTE/PCA data we built a set of source spec- 
tra in the 4 — 25 keV energy band and approximated them 
in a model that includes interstellar absorption, powerlaw 
and black-body components (wabs, powerlaw and bbodyrad 
models in the XSPEC package). Interstellar absorption in 
the direction of the source is quite low (Nn = few x 10 21 ) 
cm -2 and cannot be determined correctly with the PGA in- 
strument. It was fixed at the v alue of Nu = 3.4 x 10 21 cm" 2 
as me asured by XMM-Newton (|La Palombara fc Mereghettil 
l2007h . Dependencies of the black-body temperature kT, the 
radius of the emission surface of the thermal component 
R, the photon index V and the normalization of the non- 
thermal component A with time are shown in Fig[3] along 
with the changes of the source flux in the 2 — 20 keV en- 
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ergy band, which was calculated from the best fit spectral 
models. 

It is clearly seen that the photon index remains very 
stable during the whole period of observations and its value 
r ~ 2 is in a good agreement with broadband spectral mea- 
surements performed with the INTEGRAL observatory (see 
below). At the same time, the normalization of the non- 
thermal component A increases with the source brightness 
by a factor of 3 — 4. This increase is more significant in 
comparison with the extension of the emission surface of 
the thermal component from R ~ 110 m to R ~ 150 m that 
should lead to the increase in the relative contribution of the 
non-thermal component to the total source spectrum with 
the source brightening. But quantitatively this increase is 
practically negligible due to the simultaneous growth of the 
black-body temperature kT from ~ 1.6 to ~ 1.8 keV (Fig[2J), 
that in turn leads to the growth of the thermal component 
flux. 

It is important to note that measured values of the ther- 
mal component parameters (black-b ody temperature an d 
radius) agree with previous results of ICoburn et all (|200ll) , 
but are different from the ones obta ined with XMM-Newton 
(|La Palombara fc Meregh ctti 2007j) - the temperature is 
slightly higher, the radius is lower by a factor of two. This 
difference is due to the different energy bands of the ob- 
servatories: 0.3 - 10 keV for XMM-Newton and 4-25 keV 
for RXTE. Naturally, the capabilities of the XMM-Newton 
observatory are much more suitable for a correct determina- 
tion of quantitative characteristics of the thermal emission 
with temperatures of kT ~ (1-2) keV. Therefore the RXTE 
data reflect only the qualitative behavior (relative changes) 
of its parameters as the source flux increases. 

In addition to the long-term changes, the source X-ray 
flux demonstrates variability on the time scale of several 
days, which became especially significant (like flares) with 
the growth of the s ource average flux (Fig[5] b ottom panel, 
see also Fig.2 from iDelgado-Marti et al.ll200ll ). Such flares 
are accompanied with the variability of the source spectral 
parameters, mainly its thermal component. Note that flares 
with changes of the spectrum hardness were observed ear- 
lier during the high lumin osity states in 1970s (|White et al.l 
Il976l : iFrontera et allll979h . 



3.4 Hard X-ray emission and cyclotron line 

X Persei was observed by the INTEGRAL observatory with 
a very long exposure in 2008-2010. Taking into account that 
1) the source intensity didn't change significantly during 
these observations (Fig[T]) and that 2) the hard non-thermal 
spectral component is very stable even when the source flux 
is variable (see previous section), we could use all these 
data to reconstruct the broadband spectrum of X Persei. In 
the resulting spectrum, the signal from the source is signifi- 
cantly detected up to 160 keV (Fig(4£i) for the first time. To 
date, hard X-ray emission at energies > 100 keV has been 
detected from onl y one other low lumino sity X-ray pulsar 
RXJ0440.9+4431 (|Tsygankov et alj|2012h . 

As a first step the source spectrum was approximated by 
the powerlaw model with a high-energy cutoff that is typical 
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Figure 4. Broadband energy spectrum of X Persei obtained with 
the JEM-X (crosses) and IBIS (dark circles) telescopes of the 
INTEGRAL observatory. Solid line represents the best-fit model 
with the cyclotron absorption line in the form of cyclabs (a). 
The residuals to the fit in units of a without (V) and with (c) a 
cyclotron absorption line (see text for details). 

for X-ray pulsars \ White et al.l 1 1983T F1. However this model 
gives an unacceptable value of \ 2 = 275 for 45 d.o.f. Devi- 
ations of the measured spectrum from the model (Fig(4p) 
demonstrate a wide absorption feature near 30keV, which 
can be interpreted as a cyclotro n resonance absorpti on line 
and which was found earlier by ICoburn et al.l <|200lh based 
on the RXTE data. 

An additional component in the form of gabs (XSPEC): 

exp ^-^=pl— J exp ^ ^ E 2 <7^ CyC ^ ))' i m P roves the fit sig- 
nificantly to x 2 = 58.2 for 42 d.o.f. The resulting parame- 
ters of the cyclotron ine - the line energy E cyc i — 29.0 ± 1.0 
keV and the line width o = 9.9 ± 1.0 keV - are in a good 
agreement with the results o f lCoburn et all (|200ll ). who also 
used the gabs model. Using another suitable model cyclabs 

(XSPEC): exp I ^_ E y JZ_lU ) for the description of 

the cyclotron absorption line improves the fit even more to 
X 2 = 43.7 for 42 d.o.f. (Fig(4j:); however the line energy is 
slightly lower in this fit (~ 23.5 keV), typical for a cyclotron 
absorption line that is approximated by gabs or cyclabs mod- 
els. 

Thus, the parameters of the best-fit model (with the 
cyclabs component), describing the broadband spectrum of 
X Persei in the 4 — 200 keV energy range are: photon in- 
dex F = 2.05±g;g|, cutoff energy E cut = 69^3 keV, folding 

4 There is no necessity to include into the model an interstel- 
lar absorption and a thermal component with the temperature of 
kT ~ 1.5 keV as the sensitivity of the JEM-X telescope is much 
lower in a comparison with PCA/RXTE. Moreover, due to prob- 
lems with the JEM-X response matrix energy channels around 
6-7 keV were excluded from the analysis according to recommen- 
dations of the ISDC helpdesk (private communication) 
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Figure 5. Optical light curve of XPersei in 1996-2011 in V- 
band (top panel) along with the ASM/RXTE light curve (bottom 
panel). 



energy Ef u = 34.012° keV, line energy E, 



- 93 ^+ 1 ' 5 

cyci — ^'iJ_i 4 

keV, line width a = 13.41^ keV, line depth r = 0.56l°;oi 
Based on the measured value of the cyclotron line energy we 
can estimate the magnetic field on the neutron star surface 
as 



Bns = 



E 



2GM_ns 



11.6 



(2.4 - 2.9) x 10 12 G, 



depending on the model used (gabs or cyclabs) and us- 
ing M ns = 1.4Mq and Rns = 15 km l|SuIeimanov et all 
l201ll ) for the neutron star radius and mass estimates, re- 
spectively. We expect that at the low luminosity of ~ 10 35 
erg s _1 the accreted matter is falling down to the neutron 
star surface, where most of the energy is r eleased and the cy- 
clotro n absorption line is formed (see, e.g. iBasko fc Sunvaevl 
1 19761 ). Finally, it is important to note that the cyclotron 
line in XPersei was found to be signif icantly broader than 
typically observed in X-ray pulsars (|Coburn et al.l |2002| ; 
iFilippova et alj|2005l ). The unusual broadness may be ex- 
plained by insufficient knowledge of spectral continuum in 
the broad energy band. In particular, approximation of the 
spectrum with other components can lead to the artificial 
deficit of photons in the regio n where the spectral compo- 
nents overlap (see, e.g., Fig.8 in lDi Salvo et al]|l998l . or cur- 
rent paper of Doroshenko et al. 2012). This region roughly 
corresponds to the energy of the cyclotron line that can lead 
to the distortion of its parameters, particularly the width 
and depth of the line. 



3.5 Correlation of X-ray and optical light curves 

Since its discovery XPersei has been regularly observed in 
different wavebands. Particularly, it was found that long- 
term fluctuations in the optical brightness (V ~ 6.7—6.1) are 
accompanied by the cyclic variability in the peak ratio of H tt 
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curves. 



6. Correlation function between optical and X-ray light 



and He I lines w i th pe riods varying between ~ 1 — 10 years 
(see IClark et al.l l200ll . and references therein) . The optical 
variability of the source is usually interpreted as a change of 
the mass content and distribution in the circumstellar disk 
due to the presence of one armed density wave. The lack of 
correlation between the strength of the line and continuum 
emission demonstrates that the disk variability involves a 
redistribution of material withi n the disk, leading to changes 
in the radial density gradient (|Clark et al.ll200ll ). 

Previous X-ray observations showed a number of cases 
when a rise of the optical brightness of the source wasn't 
followed by a corresponding rise of the X-ray flux. This 
was explained by the truncation of the Be star decretion 
disk by the tidally-driven eccentric instability (see, e.g., 
lOkazaki fc Negueru cla 2001). Truncation at the 3:1 resonant 
radius produces a wide gap between the disc outer radius 
(0A7a x , a x is an orbital separation ) and the periastron dis- 
tance (0.89 a x ) ()Clark et al.ll200l"h . This wide gap depress 
the accretion even at periastron, nat urally leading to a low 
X-ray luminosity. In addition to this, IClark et al.l (|200lh re- 
port an episode of complete disc loss during the period 1988 
May - 1989 June. The time needed for the outer radius disk 
to reach the truncation radius was estimated to be about 
~ 12(a/0.2) _1 yr after the disc forma tion begins (a is the 
Shak ura-Sunyaev viscosity parameter, IShakura fc Sunvaevl 
Il976h . Thus one can hardly expect any X-ray - optical cor- 
relation before the year 2000. 

Regular monitoring of the system by ASM/RXTE since 
1996 gives us a unique chance to test this hypothesis and to 
investigate the X-ray - optical dependance. Fig[5] shows the 
evolution of the system magnitude in 1^-band (upper panel) 
along with the X-ray evolution (bottom panel). It is clearly 
seen that the large variability of the source magnitude has 
only a weak reflection in the X-ray flux before 2002, when 
the rise of the optical flux is followed by a major outburst 
in X-rays. Given the discussion above, this rise in the X-ray 
flux can be interpreted as the interaction of the neutron star 
with the fully recovered disk. The maximum rise of the X- 
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ray flux in 2003 is followed by a rapid decrease to the level 
still two times higher than the pre-outburst level, while the 
optical flux stays at a rather high level until a small drop in 
20f0 and a subsequent rise in 2011. The drop of the X-ray 
flux in 2004 is probably related to the disruption of the outer 
layers of the disk, and the second rise of the X-ray flux is 
in good agreem ent with the estim ated 5.7 years for the disk 
growth peiord (|CIark et alj|200ll ). 

In order to test qualitatively the correlation between 
optical and X-ray light curves with a possible non-zero 
time lag, we performed an interpolate d correlation analy- 
sis usi ng th e improvements proposed by lGaskell fc Peterson! 
l | 1987ft and IWhite fc Peterson! | 19941 ). In order to increase 
the statistic we have binned ASM/RXTE data into 10 day 
bins. The resulting Interpolated Cross-Correlation Function 
(ICF) is shown in Fig(5] These data show a tentative cor- 
relation between X-rays and optical flux with a zero time 
lag (with a correlation coefficient R — 0.53 ± 0.01). The 
pronounced peaks at time lags of ±1, ±2, 3, 4 years are most 
probably artificial ones due to the internal ASM/RXTE one 
year periodicity relate d to the motion of Earth around the 
Sun (|Wen et alj 12006ft and also the sampling of the opti- 
cal observations after 2003 (see Fig(6|. These facts can also 
lead to the widening of the ICF, thus despite the quite large 
correlation coefficient we prefer to call this correlation a ten- 
tative one. 



4 CONCLUSIONS 

We presented here the comprehensive temporal and spec- 
tral analysis of the X-ray pulsar X Persei using data of the 
RXTE, INTEGRAL and XMM-Newton observatories. The 
main results can be summarized as: 

- a strong outburst activity was detected from the 
source in soft and hard X-rays during 2001-2011. During this 
period the source luminosity twice reached its maximum for 
the last ~ 30 years, L x ^ 1.2 x 10 35 erg s -1 ; 

- this activity was accompanied by drastic changes of 
the pulse period: a ~ 35-years period of the spin-down 
changed to spin-up with an average rate of P/P ~ —3.6 x 
10 -4 yr _1 , several times higher than previous rates of spin- 
up and spin-down; moreover, this acceleration of the rota- 
tion of the neutron star is apparently not in steady state, but 
depends on the source flux; the observed behaviour of the 
pulse period and its quantitative properties are naturally ex- 
plained in terms of the current theory of the quasi-spherical 
accre tion in low-luminosity X-ray pulsars jShakura et al.l 
[2013); 

- at the same time, no significant variability of the spec- 
tral parameters with the source flux was found during the 
rising part of the first outburst in 2001-2003; the source 
spectrum in the 4 — 25 keV energy band can be well de- 
scribed by a combination of a thermal component (with a 
temperature of kT ~ (1.6— 1.8) keV) and non-thermal com- 
ponent (powerlaw with the photon index of F ~ 2); 

- using deep observations of the INTEGRAL observa- 
tory we reconstructed a broadband spectrum of X Persei in 
the 4 — 200 keV energy range and, for the first time, signif- 
icantly detected hard X-ray emission from the source up to 
~ 160 keV; 

- a significant detection of the cyclotron absorption line 



in the source spectrum allowed us to estimate the magnetic 
field strength on the neutron star surface as Bns — (2.4 — 
2.9) x 10 12 G depending on the model used (gabs or cyclabs); 
at the same time an insufficient knowledge of the spectrum 
continuum in the broad energy band can lead to distortion 
of the cyclotron line parameters (particularly its width and 
depth may be affected); 

- the observed variations of the X-ray flux are in good 
agreement w ith predictions of the viscous decretion disk of 
the Be star (|Clark et al.ll200lT ). a quantitative study of the 
X-ray - optical cross-correlation shows only a tentative cor- 
relation with a zero time lag. 
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